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The ionic adsorption around a weakly charged spherical colloid, immersed in size-asymmetric 1:1 
and 2:2 salts, is studied. We use the primitive model of an electrolyte to perform Monte Carlo 
simulations as well as theoretical calculations by means of the hypernetted chain/mean spherical 
approximation (HNC/MSA) and the unequal-radius modified Gouy-Chapman (URMGC) integral 
equations. Structural quantities such as the radial distribution functions, the integrated charge, and 
the mean electrostatic potential are reported. Our Monte Carlo "experiments" evidence that near 
the point of zero charge the smallest ionic species is preferentially adsorbed onto the macroparticle, 
independently of the sign of the charge carried by this tiniest electrolytic component, giving rise to 
the appearance of the phenomena of charge reversal and overcharging. Accordingly, charge reversal 
is observed when the macroion is slightly charged and the colons are larger than the counterions. 
In the opposite situation, i.e. if the counterions are larger than the colons, overcharging occurs 
(a feature originally predicted via integral equations in J. Phys. Chem. B 108, 7286 (2004), for 
the planar electrical double layer). In other words, in this paper we present the first simulational 
data on overcharging, showing that this novel effect surges, close to the point of zero charge, as a 
consequence of the ionic size asymmetry. We also find that the HNC/MSA theory captures well the 
charge reversal (CR) and overcharging (OC) phenomena exhibited by the computer "experiments", 
especially as the macroion's charge increases. On the contrary, even if URMGC also displays CR 
and OC, its predictions do not compare favorably with the Monte Carlo results. Further, it is seen 
that the inclusion of hard-core correlations in Monte Carlo and HNC/MSA leads to spatial regions 
near the macroion's surface in which the integrated charge and/or the mean electrostatic potential 
can decrease when the colloidal charge is augmented and vice versa. These observations aware about 
the interpretation of electrophoretic mobility measurements using the standard Poisson-Boltzmann 
approximation beyond its validity region. 

PACS numbers: 61.20.-p, 61.20.Gy, 61.20.Ja, 61.20.Qg. 



I. INTRODUCTION 

It is widely known in physical-chemistry that a sur- 
face in contact with an electrolyte solution usually be- 
comes charged and, thus, that the ions around the in- 
terface acquire a diffuse structure commonly denoted as 
the electrical double layer (EDL). One of the most suc- 
cessful early theories used to describe these systems in 
the dilute and/or weak electrostatic regimes is the clas- 
sic Poisson-Boltzmann (PB) treatment, which is based 
on a point-ions representation of an electrolyte [J, d, Q ■ 
Under this approach, it is an accepted fact that the coun- 
terions of a binary electrolyte are mostly adsorbed to the 
electrode when its surface charge density is increased; the 
colons, on the other hand, are pushed away from the re- 
gion close to the charged surface [1, As a result, in 
the PB picture, the main role of the EDL is to neutral- 
ize monotonically the surface charge as a function of the 
distance to the macroparticle, leading to a screened in- 
teraction between charged colloids in solution [5]. How- 
ever, starting from the middle of the past century, an 
appreciable number of experimental instances has been 
detected in which the effective charge of a colloid seems 
to reverse its sign (recent electrophoresis papers can be 



consulted in Refs. d, HI!])- This singularity, due to an 
excess in the counterion's compensation of the bare sur- 
face charge, is known as charge reversal (CR), and the 
too simplified PB formalism can not describe it since its 
occurrence involves strong ion-size correlations. It must 
be noted that, despite a lot of experimental research, 
it was only until recent times that direct measurements 
of CR were performed by Besteman et al. [1, [13, [HI, 
who employed atomic force microscopy techniques and a 
new electrophoresis capillarity apparatus to achieve lab- 
oratory conditions that were very difficult, or even im- 
possible, to reach in traditional static and electrokinetic 
experiments. 

Through the years, the experimental advances have 
prompted the development of theoretical explanations 
for CR. In this way, the strong correlated liquid for- 
malism [12| was conceived as an effort to comprehend 
and predict this kind of behavior. For macroparticles 
dissolved in trivalent electrolytes this scheme yields in- 
formation consistent with observational data, but breaks 
down when applied to divalent solutions 11]. Alterna- 
tively, the theory of integral equations for liquids has 
proved to be a robust and reliable approach in this area 
of research [U, [H, [H, [H [H, [13, El [H HI- For in- 
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stance, it has been the main route to demonstrate that 
the ionic excluded volume constitutes, by itself, a physi- 
cal mechanism with the ability to induce charge reversal 
(2]| . although this does not mean, naturally, that the CR 
observed experimentally under very diverse conditions is 
due purely to such steric effects (an interesting discus- 
sion about t he p hysical and/or chemical origins of CR is 
available in [23|). 

Another fascinating phenomenon that has been theo- 
retically predicted to occur in the EDL is overcharging 
(OC). OC appears when the colons arc adsorbed to the 
electrified interface, despite the coulombic repulsion, in- 
creasing the native surface charge. This anomaly was first 
observed and defined for a mixture of macroions and a 
size-symmetric electrolyte in contact with a charged wall 
[2^. Notably, in such article the goal was to establish 
that the overcharging was prompted by the macroions, 
whereas here we will show that OC can be present even 
in a simpler system (without a wall) if ionic size asymme- 
try comes into play. In this context, it is worth to recall 
that the most basic way to incorporate consistently ionic 
size contributions into the EDL is by using the so-called 
primitive model (PM) of an electrolyte, wherein the ions 
are taken as hard spheres with point charges embedded at 
their centers. A large amount of work has been reported 
for the EDL using the PM in the special case of equal- 
sized ions (i.e. the restricted primitive model (RPM)), 
which implies a considerable simplification at the level of 



the model and calculations 




M 1221, |28|, |29|, |30|, |31|, |3a, 1331 . Notwithstanding, it is 

more realistic to expect EDL systems with different ionic 
sizes in their electrolytic species (due, for instance, to a 
distinct degree of hydration of the ions) and, thus, the 
PM should be preferred for a more faithful representation 
of the EDL. 

Along these lines, the issue of the practical relevance of 
ionic size asymmetry in the physical-chemistry of surfaces 
has been recently revivified in an experimental paper 
on the electrokinetics of uncharged colloids by Dukhin 
et al. [sj]. In such investigation, the authors revisited 
the idea that "...a double layer might in fact exist, even 
when there is no electric surface charge at all (on the 
colloid), solely because of the difference in cation and 
anion concentrations within the interfacial water layer..." 
and provided a measurement technique and experimental 
data supporting the existence of the so-called zero surface 
charge double layer, a concept introduced in a theoretical 
model by Dukhin himself and coworkers more than two 
decades ago [s^ [s^ [s^ • The relevant fact to our discus- 
sion is that, as proposed by Dukhin and other authors 
[3^ . [39I |40| . l4ll . [43 , |43| | , such a difference in cation and 
anion concentrations, and the concomitant charge sepa- 
ration in the proximity of an uncharged colloid, can be 
attributed to the difference in the distances of closest ap- 
proach of counterions and coions to the surface (as an al- 
ternative, or in addition, to the "chemical" phenomenon 
of specific adsorption). Thence, the work by Dukhin et 
al. emphasizes the importance of the ionic size asymme- 



try effect in relevant phenomena occurring in real EDL 
systems, like the binding of simple inorganic electrolyte 
ions on colloidal substrates, and electrokinetics. 

It is interesting that in the literature there are rel- 
atively few publications considering ionic size asymme- 
try (see, for example, 0, 0, [H, gl, \M, H \M, M, M, M, 
m, [5JI), most of them dedicated to the planar geometry 
and only two reports for the spherical case . In our 

opinion, this apparent lack of interest in size-asymmetric 
electrolytes might arise from the long-standing belief in 
the dominance of counterions in the EDL; a fact foreseen 
and corroborated in a couple of pioneering papers on the 
modified Gouy-Chapman (MGC) theory for planar inter- 
faces [13, E^. To be more explicit, the precise meaning 
of the dominance of counterions is that "...at large po- 
tentials or charge densities, the coions are excluded from 
the vicinity of the electrode. Consequently, the counte- 
rions dominate and the double layer properties approach 
those of a symmetric electrolyte whose charge and diam- 
eter are equal to those of the counterions..." [45]. No- 
ticeably, and even if this behavior was originally enunci- 
ated and verified exclusively at the MGC level, during the 
past years most of the modern studies of the EDL, that 
use theoretical approaches which surpass the punctual- 
ions PB treatment, have subscribed (or assumed with- 
out a rigorous proof) such counterion predominance in 
the primitive model [H, [H, [13, H Hi [H, [H, [13, a 
situation that has resulted in the mentioned scant at- 
tention to size-asymmetric EDLs. In contrast, very re- 
cently, some of the present authors have published an 
integral equation and simulational analysis of the size- 
asymmetric spherical EDL Q where it has been evi- 
denced that, for highly charged surfaces, counterions do 
not always dominate, i.e. that coions really matter in the 
double layer. At this point, it should be also noted that 
the establishment of the counterion dominance by Val- 
leau, Bhuiyan and coworkers [3,14^ was not really based 
on the plain Poisson-Boltzmann equation, or equivalently 
on a model of a genuine punctual electrolyte, since for 
these authors the ion-ion potential corresponds to that 
between point charges whereas for the ion-wall interac- 
tion a closest approach distance (hard-core or Stern cor- 
rection) is added. In fact, the unequal radius modified 
Gouy-Chapman (URMGC) approach of Refs. [H, El 
represents not only the inclusion, at the lowest order, of 
excluded volume contributions into the Gouy-Chapman 
theory (via the Stern modification), but also the first 
attempt to take into account the ionic size asymmetry 
through the use of different distances of closest approach 
for counterions and coions. From all the above discussed, 
the new integral equations results for a colloid in con- 
tact with a size-asymmetric PM electrolyte [H, [sil] (in 
which hard-core contributions are consistently embodied 
in both the ion-surface and ion-ion interactions) imply 
that an incomplete consideration of excluded volume and 
size asymmetry contributions, like that of the URMGC 
theory, can lead to an inaccurate description of the dou- 
ble layer at high surface charges. 
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On the other hand, apart from the dominance of coun- 
terions at large electric fields, URMGC has predicted 
other notable phenomena in the EDL, this time for low 
charged surfaces, such as the appearance of a potential 
of zero charge, oscillations in the ionic density and mean 
electrostatic potential functions, and indeed charge rever- 
sal and overcharging [sl, . Lately, several theoretical 
and simulational investigations [51, 53, 60] have pondered 
the planar EDL for a PM electrolyte in the low-charge 
regime and have corroborated the existence of the poten- 
tial of zero charge and the non-monotonic behavior of the 
ionic radial distribution functions (RDFs) and potential 
profiles, previously seen in the "semi-punctual" URMGC. 
Nevertheless, in the same reports it was also found that 
only those theoretical formalisms (e.g. the MPB5 and 
density functional theories) that fully include the hard- 
core and ionic size asymmetry effects succeeded in de- 
scribing quantitatively the EDL near the point of zero 
charge [52, 3Si, contrasting with URMGC that showed 
solely a limited success for 1:1 salts (U- In addition, 
and with respect to CR and OC, in a 2006 paper [5^ Yu 
et al. noticed, for the first time, the appearance of these 
"anomalies" in an URMGC treatment of an electrolyte 
in a charged slit; a certainly intriguing fact given that, in 
the past, CR and OC had been observed only in theoreti- 
cal analysis of primitive model EDLs [2l|, [23|. Therefore, 
at present, an exhaustive simulation study that confirms 
and characterizes the phenomena of charge reversal and 
overcharging in slightly charged PM EDLs, as well as an 
application of reliable theories in order to explain these 
striking features, is still lacking in the literature. Pre- 
cisely, the main objectives of this communication are, in 
the first place, to supply fresh and comprehensive Monte 
Carlo (MC) data about the potential of zero charge, 
charge reversal, overcharging and the behavior of diverse 
structural properties of a low charged primitive model 
EDL in spherical geometry and, secondly, to present the 
comparison of such simulational information with the 
corresponding theoretical results of the HNC/MSA and 
URMGC integral equations, trying to assess the conse- 
quences of a consistent treatment of the excluded volume 
and ionic size asymmetry contributions in the spherical 
electrical double layer. 

To investigate the static properties of the size- 
asymmetric spherical EDL in the weakly charged regime, 
with special focus on the CR and OC phenomena, we 
have produced simulation and HNC/MSA integral equa- 
tion results for the radial distribution functions (RDFs) 
of 1:1 and 2:2 primitive model electrolytes bathing a 
spherical colloid under diverse conditions. From the 
RDFs we extract the mean electrostatic potential and 
the charge profiles to identify the presence of charge re- 
versal and overcharging and to examine their dependence 
on the ionic size asymmetry and other system parame- 
ters. All this structural information is contrasted with 
that corresponding to the semi-punctual URMGC the- 
ory. As it will be illustrated below, most of the computer 
"experiments" data are better paralleled by HNC/MSA 



than by URMGC. Additionally, we would like to note 
that, to the best of our knowledge, this work reports 
the first confirmation, via simulations, of charge rever- 
sal and overcharging in the EDL near the point of zero 
charge. This paper is organized as follows. In Section 
II we present the model system and the theoretical ap- 
proaches. Section III contains the details of the numer- 
ical solution of the integral equations and of the Monte 
Carlo simulations. Section IV is devoted to the results 
and their discussion, and we close with a summary of 
relevant findings and conclusions in Section V. 

II. MODEL SYSTEM AND THEORY 

Our basic representation of the spherical electrical dou- 
ble layer (SEDL) is constituted by a rigid, charged spher- 
ical colloid of diameter Dq and uniform surface charge 
density cto, surrounded by a continuum solvent of dielec- 
tric constant e. The macroion is in contact with two ionic 
species, which in the primitive model are treated as hard 
spheres of diameters Di [i = 1,2) with embedded point 
charges, g^, at their centers. Note that qiq2 < 0. The 
interaction potential between the particles in this model 
(i.e. macroion and electrolytic ions) is then given by 

[ gi(7j/(er-), for r>Dij, ^' 

where the subscripts i and j run from to 2, r is the 
center-to-center distance between two particles of types 
i and j, Dij — {Di + Dj)/2, qi — Zie is the charge of the 
species i with valence z^, e is the protonic charge, and, 
for the spherical colloid, q^ = z^^e = 47r(Do/2)^cro. The 
system as a whole is electroneutral, i.e. -^iPj = ^^ 

where pj is the bulk number density of the electrolytic 
species j. 

The Ornstein-Zernike equation for a multicomponent 
mixture of S species is [6l| 



h^3{r) = c,j(r) + ^ Pfe / h,k{t)ckj{\ r- t\)dV, (2) 

such as Pi is the bulk number density of each one of the 
species in the system, r = | r | and hij (r) are, respectively, 
the distance and the total correlation function between 
two particles of types i and j, gijir) = hij{r) -I- 1 is the 
radial distribution function, Cij (r) is the direct correla- 
tion function, and i = 1 1 1 being the distance between 
two particles of types i and k. This group of equations 
requires a second relation (or closure) between the total 
and direct correlation functions. For charged systems, 
the hypernetted chain (HNC) and the mean spherical 
approximation (MSA) closures are widely used [6ll . [63 | . 
The HNC and MSA relations, for r > Dim, are given as: 
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cini{r) = -f3Uim{r) + hira[r) - \n[hhn{r) + 1], (3) 
for HNC, and 

Ci^[r)^~(iUi^{r), (4) 

for MSA, where (3 = l/{kBT) is the inverse of the thermal 
energy. These expressions are complemented by the exact 
condition hi,n{r) = — 1, for r < Dim- 
Let us consider 5 = 3 and that the species corre- 
sponds to macroions (thereinafter denoted equivalently 
by M) at infinite dilution, whereas the indices 1 and 2 
are associated to a binary electrolyte. Then Eq. ([2]) for 
the components M{= 0) and j can be written as: 



hM]{r) ^ CMj{r) +'^Pk hMk{t)ckj{\r-t\)dV, (5) 
k=i 

Note that Eqs. ([5]) are a complete set of integral equa- 
tions for the SEDL. When Eq. ([3]) is employed in Eqs. ([5|) 
solely for CMj{r), and the Ckj{\r — t\) are approximated 
by the MSA analytical expressions for a bulk electrolyte 
[si il the HNC/MSA integral equations are ob- 
tained. A detailed account of this HNC/MSA formal- 
ism can be consulted elsewhere ^ and will not be re- 
peated here. Our election of the HNC/MSA theory to 
perform the present study is based on the fact that, for 
many years, it has been used successfully to investigate 
the RPM electrical double layer in many geometries (e.g . 
planar, cyhndrical, and spherical) [sl. IT4. [l6l. ITtI. l66l 

The integral version of the URMGC theory in spherical 
geometry is easily deduced from the HNC/MSA formu- 
lation if Ckj{\ r — 1 1) — ^PqkQj / {e\ f ~ t\) is inserted in 
Eqs. (O, instead of the inter- ionic MSA direct correlation 
functions. It must be stressed that, in the original pa- 
pers m,!!^, as well as in all the posterior treatments [Slj, 
URMGC has been solved, in differential form, strictly for 
planar interfaces, which means that the present study ex- 
tends the classic URMGC planar theory to the spherical 
instance, continuing along the lines of our previous UR- 
MGC (and HNC/MSA) account of the SEDL [5;]. 

Once the gMjii") are available, either from a theory 
(e.g. HNC/MSA or URMGC) or from simulation, it 
is then possible to calculate various relevant functions, 
namely the local electrolyte charge density (LECD), 

2 

P* (^) = 51 ZjPj9Mj (r)47rr^ , (6) 

i=i 

the total integrated charge (IC), 

P{r)^ZM+ I P*{t)dt, (7) 
Jo 



and the mean electrostatic potential (MEP), 

m--J^^dt. (8) 

These quantities are fundamental in our analysis of the 
properties of the SEDL. The LECD is a linear density 
that gives us information about the electrolytic charge 
(in units of e) inside a spherical shell of thickness dr lo- 
cated at a distance r from the center of the macroion. 
Besides, the integral over all the space of —p*{r) results 
in the valence of the macroion, zm, as required by the 
electroneutrality condition. On the other hand, the IC is 
a measure of the net charge (in units of e) enclosed in a 
sphere of radius r centered in the macroion. It equals zm 
for Dq/2 <r < (Do + Di)/2, if Di < D2, and goes to 
zero as r — !■ 00, again due to the electroneutrality restric- 
tion. The IC also allows to compute the amount of charge 
adsorbed onto the macroparticle, i.e. the accumulated 
charge within the Helmholtz planes (see below for the 
definition of the Helmholtz planes), and, moreover, it de- 
tects charge reversal when P{r)z},[ < 0, and overcharging 
if P{r)zM > and \P{r)\ > \zm\- In addition, the MEP 
around the macroion is a central magnitude in colloid 
science because it determines, for instance, the regimes 
of stability/flocculation or the migration of macroparti- 
cles in a colloidal suspension [6^. As a matter of fact, 
the MEP at certain distance near the macroion's surface 
is usually identified with the well-known clcctrokinetic 
potential at the slipping plane (or zeta potential, 
The C potential is experimentally measurable in systems 
that display electrokinetic effects such as electrophore- 
sis, electro-osmosis and streaming currents, and it has an 
ample use in the physico-chemical characterization, sep- 
aration and/or fabrication of colloidal materials 0, [6^. 

In particular, we will be interested in the IC and the 
MEP in the neighborhood of the Stern layer. More es- 
pecifically, the Stern layer is the free-of-ions space next 
to a macroion that ends at the Helmholtz plane. The 
Helmholtz plane (or, more properly, the Helmholtz sur- 
face) is the geometrical place associated to the colloid- 
ion closest approach distance ,68.] . In size-symmetric 
electrolytes only one Helmholtz plane can be identified. 
In our model, however, the size asymmetry between the 
ions allows us to define an inner Helmholtz "plane" (IHP) 
and an outer Helmholtz "plane" (OHP) (note the conven- 
tional usage of the word "plane" ). The IHP is specified by 
the closest approach distance of the smallest ionic com- 
ponent to the colloidal surface (i.e., by {Dq + Di)/2, if 
Di < D2), whereas the OHP is determined by the cor- 
responding distance of closest approach for the largest 
species (i.e., by (Dq + D2)/2, if Di < 02). Therefore, 
for the primitive model EDL, if Di < D2, the Stern 
layer is the region where Do/2 < r < {Dq + Di)/2, and 
the Helmholtz zone corresponds to {Dq + Di)/2 < r < 
{Do + D2)/2. Obviously, for size-symmetric salts the IHP 
and the OHP coincide and the standard notions of the 
Helmholtz plane and Stern layer are recovered. In the 
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general PM case, and provided that Di < D2, when the 
MEP is evaluated at r = {Dq + Di)/2, Eq. (O gives 
the MEP at the IHP, which we denote as tpiHP- On the 
other hand, if Eq. ^ is calculated at r = {Dq + D2)/2, 
the MEP at the OHP, ipoHP, is obtained. 



III. NUMERICAL DETAILS 

A size-asymmetric 1:1 or 2:2 electrolyte with a ratio 
between ionic diameters D-/D+ = 2, bathing a charged 
macroparticle of diameter Dm = Dq = 20 A and va- 
lence zm, was considered in all the calculations reported. 
Specifically, the diameters of the positive and negative 
species were D^ = 4.25 A and D^ = 8.5 A, respectively. 
In other words, for definitiveness, the cations have been 
chosen as the smallest species of the binary electrolyte, 
i.e., — Di and D_ = D2, with I?i < D2. Note also 
that the sign of zm defines which one of the ionic species 
is the colon or countcrion. 

To establish the primitive and semi-punctual models 
employed in our simulational and theoretical approaches, 
let us introduce the macroion-ion contact distances, (i_A/+ 
and d-M-, given by 



dMi = < 



{Dm + D+)/2, 
{Dm + D^)/2, 



for i = +, in simulation, 
HNC/MSA and URMGC, 
for / = — , in simulation, 
HNC/MSA and URMGC. 



(9) 



It must be remembered that here the macroions, cations 
and anions correspond to the indices 0, 1 and 2, respec- 
tively. 

Complementarily, the ion- ion contact distances, + , 
d , and d+-(= d- +), are 



D+, for i = j = in simulation 

and HNC/MSA, 
D- , for i = ji = — , in simulation 

and HNC/MSA, 
(Z)+ -|- D-)/2, for i = + and j = — , in simulation 

and HNC/MSA, 
0, for any i and j, in URMGC. 

(10) 



The dielectric constant and temperature considered in all 
the cases were e = 78.5 and T = 298 K. 

The URMGC and HNC/MSA theories were numeri- 
cally solved by means of a Picard iteration scheme, which, 
in the past, has been thrivingly employed in a number 
of studies of the EDL in various geometries via integral 
equations and density functional theories [:5, 29i, 3>||. 

Simulations were performed in a cubic box with the 
usual periodic boundary conditions in the canonical en- 
semble. The ionic species satisfied the electroneutrality 
condition: zm + -^+-2+ + N^z^ + NcZc — 0, where zm 
is the valence of the macroion, N^, z+ and A^_, z_ are 



the number of ions and the valence of the positive and 
negative electrolytic species, respectively, and Nc and Zc 
are the number and the valence of the counterions that 
balance the colloidal charge. In order to accomplish con- 
sistent comparisons with the theory, the absolute value 
of the valence of such counterions was \zc\ — 1 for the 1:1 
electrolyte and \zc\ = 2 for the 2:2 salt. The macroion 
was fixed in the center of a simulation box of length L 
and, in order to avoid border effects, the extension of 
the cell was enlarged until the integrated charge showed 
clearly a plateau of zero charge far from the macroion. 
The runs were done for w 2000 particles for the mono- 
valent salt and w 1000 ions for the divalent electrolyte. 
The long-range interactions were taken into account by 
using the Ewald sum technique with conducting bound- 
ary conditions [toI Fn| . The associated damping parame- 
ter was a = b/ L and the fc-vectors employed to compute 
the reciprocal space contribution satisfied the condition 
fc| < 5. Monte Carlo runs of charged systems were per- 
formed with 5 X 10^ equilibration MC cycles and from 
6 X 10^ to 1.8 X 10^ MC cycles were practised for the pro- 
duction stage. The radial distribution functions were cal- 
culated using standard bin procedures [tO, [71| , whereas 
the integrated charge and the electrostatic potential were 
obtained by using Eqs. ([7]) and ([5]), respectively. 



IV. RESULTS AND DISCUSSION 

In what follows, our results are discussed chiefly in 
terms of the Monte Carlo (MC) simulations data. Ap- 
propriate comparisons with the HNC/MSA and URMGC 
formalisms are presented such that the accuracy of the 
theoretical predictions can be assessed. 



A. Monovalent size-asymmetric electrolytes 

The structure of the EDL is the result of an cntropic 
and energetic competition. In the proximities of the point 
of zero charge the entropy is expected to be important, 
whereas for highly charged macroparticles the EDL is ex- 
pected to exhibit strong coulombic correlations. Thus, in 
order to understand the behavior of the ionic atmosphere 
next to a barely charged colloid in contact with a size- 
asymmetric electrolyte, we will first present a compari- 
son between the surrounding distribution of a univalent 
salt and that of a mixture of hard spheres having the 
same radii and concentration. This will illustrate how 
the structure of charged systems deviate from the neu- 
tral situation. 

Let us consider initially the ionic distribution around 
an uncharged macrosphere (zm = 0). Fig. [1] (main 
panel) displays the radial distribution functions (RDFs) 
of two systems; in one case the ionic species represent a 
1:1, 1 M salt, and in the other the ions are uncharged, 
forming a pure hard-sphere assembly. We note here that, 
in the former system, the EDL is the result of both 
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the entropic and energetic contributions; in contrast, the 
structure of the latter instance is only driven by entropy 
because all the interactions are of the excluded volume 
type. Nevertheless, since we are working in the zero col- 
loidal charge regime, we expect that in the first case 
the excluded volume interactions must play a determi- 
nant role in the resulting properties of the double layer. 
The direct comparison of the simulational and theoreti- 
cal RDFs of the charged and uncharged spheres around 
the neutral macroion is presented in the main panel of 
Fig. [T] An inspection of the Monte Carlo data shows 
that the structure of the two mixtures (of ions and of the 
hard spheres) is rather similar for r'/I?_|_ > 2.5 (where 
r' is the distance measured from the colloidal surface), 
indicating the existence of weak charge correlations at 
these distances. Differences, however, are noticeable at 
smaller distances, particularly for the tiniest species. We 
observe that the addition of charge to the ionic species 
slightly increases the concentration of the larger species 
and decreases the concentration of the smaller species, 
especially at the contact distances. As it is shown later, 
this dewetting of the surface augments with the valence 
of the ions. Such comportment of the RDFs could lead 
eventually to significant changes in the thermophysical 
properties of the charged systems since, for instance, the 
pressure depends directly on the contact values of the 
pair correlation functions (6^ . We must point out that 
these are rather concentrated systems, with an ionic vol- 
ume fraction oi (j) k, 0.217. That is the reason why the 
contact peaks of the neutral system are so high. Hence, 
at these volume fractions, the charge correlations in the 
neighborhood of the macroparticle are masked by the 
strong steric contributions. Notwithstanding, we remark 
that, for barely charged systems, the absolute values of 
the ionic RDFs close to the surface and the extent of their 
deviations with respect to the pure hard-sphere mixture 
will rule the behavior of the IC and MEP and, as it is evi- 
denced below, will be crucial to the degree of appearance 
of CR and OC. 

The adequacy of a theoretical description of the 
present EDL systems will depend on its ability to capture 
correctly the electrostatic and steric correlations close to 
the macroparticle. In this regard, the integral equations 
results portrayed in Fig. [T] illustrate that HNC/MSA 
follows closely the trends of the simulations, with quan- 
titative discrepancies near the colloid where this scheme 
overestimates the RDFs. In contrast, the URMCG data 
exhibit very different tendencies from those of the simu- 
lations. Especially noticeable is the pronounced separa- 
tion between the URMGC radial distribution function of 
small ions and those from HNC/MSA and Monte Carlo, 
as well as the very low values of the URMGC normalized 
density of cations in the zone comprised by the Hclmholtz 
planes. In fact, such exaggerated absence of small ions in 
URMGC will be of consequence in our later analysis of 
charge reversal and overcharging. Also, and contrary to 
the MC data, the URMCG theory predicts g{r)s that are 
monotonic beyond the OHP, a well documented charac- 



teristic of point-ions theories 0, d, [6^ . From the previ- 
ous discussion, it is therefore expected that the ensuing 
properties of the EDL extracted from the HNC/MSA and 
URMGC ionic profiles should present important discrep- 
ancies, with HNC/MSA excelling in the comparison with 
Monte Carlo. 

To better visualize the distribution of charge around 
the neutral macroparticle, let us focus on the function 
p* (r) . The analysis of the inset of Fig. [1] shows that 
the Monte Carlo LECD around the colloid exhibits a 
structure with alternating domains of positive and neg- 
ative local charge. Three intervals are clearly seen in 
the scale of the graph. The first domain is positive and 
encompasses the region between the Hclmholtz planes. 
A second space, where p*{r) is negative, spans from the 
OHP up to r'/D+ « 2.4, and a third interval, where 
p*{r) > 0, is located at 2.4 < r'/D+ < 3.4. In the first 
domain the local charge is positive because the larger 
ionic species (the anions) is completely excluded from this 
region. Beyond the OHP p*{r') drops sharply to a large 
negative value, due to the presence of a compact contact 
layer of anions at the OHP. The charge remains negative 
in this second region until the concentration of positive 
and negative ions equalize, namely at the first crossing 
point of the corresponding radial distribution functions. 
Clearly in the third zone the population of cations dom- 
inates, leading to the positive charge observed in the 
graph. This pattern of alternating domains continues 
as one goes away from the surface of the macroion, al- 
though they are not distinguishable in the scale of the 
figure. Again, we should note that HNC/MSA parallels 
very nicely the MC simulations, whereas URMGC pre- 
dicts a monotonic, asymptotic increase of the local charge 
outside the Hclmholtz planes. Thus, for the last theory, 
the charge outside the Hclmholtz planes, at any distance 
from the macroion, is always dominated by the anions, 
which departs evidently from the predictions of Monte 
Carlo and HNC/MSA. 

Complying with the behavior of the local charge, and 
given that ^^^^ = p*{r') (see Eq.©), the P(r') variates 
as it is observed in Fig. [2] (main panel). Interestingly, we 
note that all the simulational and theoretical IC curves 
present the adsorption of a layer of charge very close to 
the macroion's surface; such layer begins at the IHP and 
reaches its maximum value of charge at the OHP. The 
amount of charge in the Hclmholtz zone is significant 
(« 3.8e,4.1e and 2.2e for MC, HNC/MSA and URMGC, 
respectively) and, in turn, gives rise to the formation of 
a double layer beyond the OHP (see Fig. [T]). The exis- 
tence of this zero surface charge double layer (ZSC-DL) 
in our calculations is regarding since confirms the ideas 
originally proposed by Dukhin et al. f33| and, in ad- 
dition, agrees with previous studies of the ZSC-DL in 
planar geometry via Monte Carlo simulations | 5 l |. and 
the URMGC and density functional theories [4^, [53| . In 
the MC and HNC/MSA cases, after its maximum, the IC 
decreases until changing sign and, subsequently, exhibits 
a region of negative values for 1.8 < r'/D^ < 3.0. Far- 
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ther than r' /D+ = 3.0 the accumulated charge fluctuates 
around zero and, finaUy, the electroneutrahty condition 
is obtained when r' oo. Contrarily, the P{r') of UR- 
MGC goes uniformly to zero. 

The information presented so far clearly evince that, at 
the point of zero charge, a very simple 1:1 salt in contact 
with a macroparticle already displays highly nonlinear ef- 
fects such as charge adsorption (i.e. a ZSC-DL) or the re- 
version in the sign of the IC. These phenomena are due to 
both the finite ionic size and to the asymmetry between 
the electrolytic species of the model. Since HNC/MSA 
fully incorporates such conditions, it is able to reproduce 
all the characteristics observed in the simulations, even 
at a quantitative level. In contrast, URMGC does not 
embody completely the ionic size correlations, just the 
non-zero contact distances between the colloid and ions. 
This is enough to capture the adsorption of charge, but 
not the additional traits of the accumulated charge at 
intermediate and large r'. In particular, URMGC fails 
to detect the sign reversal in P{r'), predicting instead a 
monotonic neutralization of the effective charge adsorbed 
inside the Helmholtz planes. In other words, the par- 
tial inclusion of the ionic size and size asymmetry con- 
tributions in URMGC has the severe inconvenience that, 
in this semi-punctual approach, the occurrence of steric- 
related pecuHarities, such as the ZSC-DL, charge rever- 
sal and the oscillation of the RDFs and IC, is restricted 
exclusively to the region between the Helmholtz planes. 
This fact, to which we will refer to as the localization of 
the size and asymmetric effects in URMGC, will be a re- 
current issue in our posterior discussions of the MEP and 
overcharging. 

The mean electrostatic potential as a function of the 
distance to the uncharged macroparticle is plotted in the 
inset of Fig. [51 The first thing worth to be noticed is that, 
despite a zero charge on the colloid, the MC, HNC/MSA 
and URMGC mean electrostatic potentials at the IHP 
are positive. This potential of zero charge (PZC) has 
been largely recognized as a direct consequence of the 
ionic size asymmetry in the EDL since the initial pa- 
pers by Valleau and others [3, SHi and, in the mean 
time, it has received great attention in diverse simula- 
tional and theoretical accounts of the planar double layer 
[13, [53. Evidently, this PZC is due to the domi- 
nant population of cations close to the macroion. Far 
from the surface of the macroion the MEP of MC and of 
HNC/MSA tend to zero, but for intermediate distances 
they have a series of minima and maxima of alternating 
sign. In particular, the first minimum defines the region 
where the MEP reversal is stronger. Now, if the wonted 
association between the electrokinetic and mean electro- 
static potentials is invoked [3], the mentioned comport- 
ment of the MEP suggests that a macroparticle could ex- 
perience electrophoresis, even if it is uncharged. Thence, 
for small cations and depending on the precise localiza- 
tion of the slipping surface [IHl , the neutral colloid should 
move in the direction of the applied field if this surface 
is somewhere in between the Helmholtz planes, or back- 



wards if the shear boundary is around the first minimum. 
It is generally accepted that the slipping or ^-plane is very 
close to the surface of the macroions [6^ , therefore we ex- 
pect the first scenario to be more plausible. Note that 
UMRGC foresees that the macroions should fiow always 
in the direction of the external electric field. 

So far we have examined systems at the point of zero 
charge, where no electrostatic correlations between the 
macroparticles and the surrounding ions were considered. 
In spite of that, the colloid- ion cntropic contributions and 
the inter-ionic correlations led to interesting phenomena. 
By weakly charging the macroions, conspicuous effects 
such as charge reversal and overcharging now arise, as it 
is seen in the remaining of this section. Then, let us eval- 
uate two situations in which the valence of the macroion 
is zm = —4 and zm = 4 (i.e. surface charge densities 
a = ±0.051 C/m^), immersed in the same 1:1 electrolyte 
as before. Note that by virtue of zm the role of anions 
and cations as counterions or coions is interchanged. Fig. 
[3] includes the radial structure (main panel) and the local 
charge (insets) of the electrolyte around the macroion for 
the two values of zm ■ Fig. [3^ contains the case in which 
the counterions are smaller than the coions (zm = ~4), 
whereas the case of larger counterions {zm = 4) is re- 
ported in Fig. [3Jd. Fig. [3^ reveals that, compared to 
the RDFs of a neutral macroparticle in Fig. [Tl the pres- 
ence of the counterions in MC and HNC/MSA is greater 
when the surface is negatively charged, as evidenced by 
the increased contact peak, whereas the concentration 
of coions diminishes. This fact stresses the relevance of 
charge correlations (induced by zm) on the ordering of 
the ions around the macroparticle. Clearly, the RDFs 
of URMGC disagree with MC and HNC/MSA. On the 
other hand, for zm = 4 (counterions larger than coions, 
Fig. [3)d), we see a dramatic decrement in the contact 
peak of the smaller ions (coions) and an augment in that 
of the larger ions (counterions), which is adequately re- 
produced by HNC/MSA and not by URMGC. Since we 
are dealing with charged systems, size and charge corre- 
lations are into play. As a result, these two cases produce 
LECD functions that look very different (see the insets 
of Fig. [3]). In addition, both insets show that the local 
charge profiles of simulation and HNC /MSA exhibit a se- 
quence of positive and negative regions. Quantitatively, 
however, the first positive zone is much more pronounced 
for zm — —4, whereas for zm = 4 the second and third 
domains are much more important. In particular, it is 
worth to notice that, for simulation and theories, the first 
region is positive in the two cases, even if the charge of the 
macroion is positive, which indicates that, for zm > 0, 
this first layer, rather than to screen the bare charge of 
the macroion, it emphasizes such native charge, which 
is a remarkable effect induced by the smaller size of the 
coions. 

Further details can be grasped by looking at the P{r') 
curves, as they are presented in Fig. 3] For zm = —4 
(Fig. 2^) we find that the Monte Carlo IC increases 
almost linearly inside the Helmholtz planes, inverting 



8 



its original sign (i.e. experiencing charge reversal) and 
reaching a maximum of Pmax{T') ~ 1.3 at the OHP. 
After the OHP, P(r') acquires an fluctuating behavior, 
where subsequent charge inversions can be appreciated. 
The inset of the figure indicates that the mean electro- 
static potential is also oscillatory, with a maximum in- 
version inside the Helmholtz zone. The level of accuracy 
of the theories can be readily estimated from the figure 
and its inset. Particularly, we see that, differently from 
HNC/MSA, URMGC is unable to describe any rever- 
sal of the accumulated charge or the mean electrostatic 
potential throughout all the space. This represents an 
extreme manifestation of the so-called localization of the 
size and asymmetric effects in URMGC. 

The functions P(r') and ip{r') for zm = 4 are por- 
trayed in Fig. and its inset. Notably, the Monte Carlo 
P{r') increases in between the Helmholtz planes, reveal- 
ing an adsorption of charge of the same sign as that of the 
macroion, which in turn "augments" the native macroion 
charge up to a maximum value of Pmax{r') k, 6.3 at 
the OHP. This striking event is referred to as over- 
charging. Overcharging (OC) was predicted theoreti- 
cally by Jimenez-Angeles and Lozada-Cassou 23| for a 
charged wall in contact with a mixture of macroions and 
a fully-symmetric salt. To this date, this peculiarity had 
not been confirmed through simulations or experiments. 
Hence, our results provide the first simulational evidence 
of OC in a very simple model, where such "anomaly" is 
a direct consequence of the ionic size asymmetry when 
charge correlations are not too strong, i.e. close to the 
point of zero charge. In this sense, overcharging appears 
when size correlations dominate over charge correlations 
near the point of zero charge. This implies, therefore, 
that OC is also expected in other geometries, whenever 
ionic size asymmetry is present. We would like to point 
out that the results presented here are in line with our 
previous reports 0, 113], where it has been established 
that the properties of the electric double layer depend on 
the two species of a binary electrolyte, and not only on 
the counterions as it has been widely accepted. On the 
other hand, regarding the performance of the theoretical 
schemes we are working with, we realize that URMGC is 
capable to describe the overcharging, but not the oscilla- 
tions and the charge reversal observed in the MC data, 
nor the inversion of the corresponding MEP. This fact 
is a typical expression of the URMGC restraint of all 
the excluded volume and size asymmetry phenomena to 
occur exclusively in the Helmholtz zone. Contrastingly, 
HNC /MSA collates very well with the simulational data, 
describing correctly all these characteristics. 

To analyze the behavior of the integrated charge as we 
depart from the point of zero charge, in Fig. [S]we present 
P(r') of MC and HNC/MSA for two sequences in zm- 
When counterions are smaller than coions {zm < 0, Fig. 
[5^) , the maximum simulational charge reversal observed 
at the OHP decreases as the valence of the macroion be- 
comes more negative. In fact, charge reversal within the 
Helmholtz planes disappear altogether for zm ~ —7. For 



Zm < —7 the first layer of charge reversal shifts further 
away from the macroparticle. On the other hand, in the 
curves corresponding to coions smaller than counterions 
{zm > 0, Fig. [Sb) overcharging occurs, and the dif- 
ference between the peak of overcharging (at the OHP) 
and Zm (at the surface) decreases as zm increases, which 
is the result of the growing macroion-coions repulsion. 
Overcharging virtually disappears everywhere when this 
electrostatic repulsion becomes strong enough, which for 
this system corresponds to zm ~ 24. Consequently, the 
maximum of overcharging only happens at the OHP, in- 
dicating that it is a feature directly caused by the size 
asymmetry of the ions in a binary electrolyte. From Fig. 
[5] the good coincidence between the simulations and the 
HNC/MSA theory is manifest. 

A more careful inspection of the Monte Carlo IC 
graphs when zm < (see the inset of Fig. [5^) shows 
that, at least for short distances, the profiles seem to 
cross each other at specific points, defining in this way 
intervals where the integrated charge as a function of zm 
follows a regular pattern, as we explain next. Thus, for 
instance, < r' / D-^- < 1.6 (note the upward arrow in the 
figure) delimits a region where, for fixed r', the simulated 
IC decreases when zm diminishes. HNC/MSA exhibits 
the same trend for < r' /D^ < 1.55 (downward arrow). 
The next crossing point of the charge profiles determines 
another space where, for a given position, the IC now 
augments when zm decreases. For simulations this oc- 
curs in the interval 1.6 < r'/D^ < 2.9 (upward arrow), 
and in 1.55 < r'/D+ < 3.1 for HNC/MSA (downward 
arrow). Below, we will return to explore the relevance 
of these findings. Previously, we proceed to consider the 
P(r') in URMGC. 

The integrated charge of URMGC does not follow any 
of the aforesaid tendencies. In the first place, for zm < 
(Fig. [6^), such theoretical approach exhibits charge re- 
versal only for a colloidal charge of zm = — 2 and this re- 
versal is weak {Pmax{r') ~ 0.7 at the OHP), which means 
that, differently from Monte Carlo and HNC/MSA, the 
CR disappears already for zm < —2. The origin of this 
comportment of the CR can be traced back to the very 
low values of the URMGC concentration of small cations 
near a discharged surface reported in the main panel of 
Fig. [1] When the negative charge on the initially neu- 
tral colloid starts to grow (i.e. zm becomes more nega- 
tive), the small (and positive) counterions are attracted 
to the macroparticle and increase their number in the 
Helmholtz zone, then producing charge reversal for very 
low Zm (—2 < Zm < 0, for URMGC). However, due pre- 
cisely to the mentioned URMGC scarcity of small coun- 
terions at Zm = 0, the growth of the negative surface 
charge overtakes that of the charge reversal owed to the 
counterions, and rapidly suppresses it. Secondly, and in 
further contrast with MC and HNC/MSA, in the data 
of Fig. [6^ we notice that, once the reversal of charge 
in URMGC ceases to occur in the Helmholtz region (for 
Zm ^ ~2), it is never observed again at any point be- 
yond the OHP in the monotonically decreasing profiles 
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of P{r'). Evidently, this is another example of the lo- 
calization of the ionic size and size asymmetry effects in 
the semi-punctual URMGC formalism. Complementar- 
ily, for ZM > (Fig. the URMGC integrated charge 
presents overcharging within the Hclmholtz planes, which 
diminishes slowly as zm augments. Notwithstanding, the 
OC in this theory is less important than that in simula- 
tions and HNC/MSA, as can be verified from a compar- 
ison of the respective differences \P{r') — zm\- Besides, 
and as expected, the IC curves in Fig. go uniformly to 
zero after the OHP. Such abatement of the overcharging 
in URMGC for positive surface charges can also be ex- 
plained in terms of the low presence of small ions close to 
a discharged macroparticle. In this case, when the posi- 
tive colloidal charge is being incremented from the zero 
value, the shortage of small cations (colons) available to 
build-up the overcharging in the proximities of the sur- 
face lessens the magnitude of the anomalous effect. 

On the other hand, apart from the usual monotonic- 
ity of the URMGC integrated charge as a function of 
the distance, in the two sequences of curves for different 
Zm presented in Fig. [51 we observe an extra monotonic 
character of the P{r') profiles, this time with respect to 
the variation of zm- Specifically, in Figs. [6^ and [6)3 the 
IC curves associated to distinct values of zm never cross 
each other. Thus, we find that the URMGC integrated 
charge functions, for any fixed r', satisfy the condition 
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> 0. 
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Otherwise, and as it was prefigured in relation with the 
inset of Fig. the integrated charges of MC and 

HNC/MSA show regions where the alternative condition 



dan 



< 



(12) 



is satisfied. We realize that this last "anomalous" behav- 
ior is possible due to the presence of crossing points in 
the corresponding IC profiles. Physically, Eq. (fTT|) estab- 
lishes that if the charge over the macroion is increased, 
the corresponding effect in the spherical EDL is to aug- 
ment locally the IC, which is intuitively awaited. Con- 
trastingly, Eq. states that the enlargement of the 
macroion's charge promotes a local decrease in the IC. 
We refer to this behavior as a local inversion of deriva- 
tive of the integrated charge (LIDIC). It should be noted 
that such inversion of the derivative occurs only when 
Zm < 0. This same phenomenon seems to be absent for 
MC and HNC/MSA if zm > (see Fig. [5]d) and, instead, 
it is apparent that all the IC profiles meet at their first 
minimum [r' / « 2.3). 

In order to explore the behavior of the potential-charge 
relationship in the SEDL, the simulational and theoreti- 
cal mean electrostatic potentials at the Hclmholtz planes 
as functions of ctq are plotted in Fig. [71 As it is usual 
in the analysis of the V'/-f/p(o'o) and tpOHpio'o) functions. 



the colloidal charge is specified here in terms of the sur- 
face density co (= ^Afe/(7rD^y^)). The reader can easily 
pass from ao to zm if realizes that the symbols (open 
circles) in the figure correspond to the sequence of in- 
teger valences zm — —8, —6, —4, —2, 0, -1-2, -1-4, -1-6, -1-8. 
In Figs[7^ and [7)3, we observe that the ipiHp{<^o) and 
ipOHpi'^o) curves for the three approaches, namely Monte 
Carlo, HNC/MSA and URMGC, display an increasing 
monotonic behavior and positive potentials of zero charge 
(PZCs). Besides, the simulations and theories predict 
the existence of intervals of negative colloidal charges, 
near the point of zero charge, for which the potentials at 
the IHP and OHP can be positive, i.e. where the condi- 
tions TpiHpo'o < and ipoHpo'o < are accomplished. 
As it is indicated by Eq. ([5]), these attributes of 4'ihp 
and ipoHP can be inferred, of course, from the comport- 
ment of the integrated charge, or better from the function 
P(r)/r^ (which is basically the local mean electrostatic 
field around the macroparticle). To facilitate our argu- 
ment let us start with the case of the URMGC theory. In 
Fig. [2]we found that, for an uncharged colloid, P{r') > 
for all r' in the URMGC scheme. Therefore, from the def- 
inition of the MEP (Eq. O), it follows that i/)(r') must 
be positive at any distance from the neutral surface; in 
particular, the relations ipiHP > and ipoHP > hold. 
This explains the positive URMGC PZCs at the IHP and 
OHP. Moreover, the series of IC curves graphed in Fig. 
[6] show a uniform precedence with respect to the varia- 
tion of (To, or that P{r' ,(Tq) > P{r',aQ) if ctq > Cg. We 
can state this equivalently by saying that there are not 
LIDIC regions in the sequence of curves in Fig. [HI This 
necessarily implies that ip{r',ao) > ^(r' ,0'^) at any dis- 
tance from the neutral macroparticle, which elucidates 
the monotonically increasing behavior of the ipiup and 
fpOHP for URMGC noticed in Fig. [71 This same prece- 
dence of P{r',ao) with respect to ctq signifies that, at a 
fixed distance from the macroparticle, there must be a 
a^"* for which 'ilj{r',a^"*) = 0. Given that in the UR- 
MGC formalism such erg'''* happens for both the IHP and 
OHP potentials within the range plotted in the figure, we 
thus find negative values of cto, close to the point of zero 
charge, associated to positive values of ipiHP and iJjohp- 

The properties of the tpiHpio'o) and ipOHpio'a) ob- 
tained from Monte Carlo and HNC/MSA can be ana- 
lyzed using similar ideas to those just applied to URMG. 
The task is more complex, however, because P(r)/r^ 
is an alternating function for both the simulations and 
HNC/MSA, which would require a more detailed exami- 
nation of the positive and negative areas involved in the 
evaluation of Eq. ([5]) . Fortunately, such process can be 
clarified if we note that: (i) in Fig. [2] the ICs of MC and 
HNC/MSA for a neutral colloid have a very high posi- 
tive peak centered at the OHP, which contributes more 
than the adjacent minimum at r' « 2.4 to the integral 
of P(r)/r^ when ipmp and ipoHP are being calculated, 
and (ii) that precisely these first maximum and minimum 
in P(r') dominate the integrals associated to ipiHP and 
ipOHP (V'/ffP-integral and '0OffP-integral, respectively). 
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From (i) and (ii) the positive sign of the Monte Carlo 
and HNC/MSA PZCs then arises. The recognition of 
the determinant role that the first extrema in the ICs 
have for the evaluation of ipiHP and ipoHP is very impor- 
tant since it provides the key to understand the complete 
shape of the ipiHpijJo) and tpoHpi'yo) coming from MC 
and HNC/MSA. The rational proceeds as follows. Dif- 
ferently from URMGC, in Figs. [5^ and[SjD the sequences 
of ICs of simulation and HNC/MSA for distinct colloidal 
charges are non-monotonic with respect to r' and dis- 
play LIDIC regions (Fig. [5^). Notwithstanding, in Fig. 
[5}d it can be seen that, for ctq > 0, the height of the over- 
charging peak at the OHP increases when ctq augments, 
or else that the fact (i) of the previous paragraph occurs 
in a more emphatic way as ctq heightens. Consequently, 
such enhancement of (i), combined with the prevailing 
fact (ii), results in the increasing behavior of ipiHp{'^o) 
for positive cro observed in simulations and HNC/MSA. 
Otherwise, for ctq < 0, it is better to consider the OHP 
as a reference point to split the contributions of P{r') to 
the V-'/fff -integral. In this form, it can be seen that, for 
Ctq ^ —0.051 C/ni^ (zji/ < —4), the area comprised by 
that section of the IC before the OHP becomes more neg- 
ative, and more important than the area enclosed by the 
IC beyond the OHP, when ctq decreases. Thus, such argu- 
ment explicates the increasing comportment of 4'ihp{o'o) 
for negative ctq, as well as the existence in Monte Carlo 
and HNC /MSA of an interval of negative surface charges 
for which ipiHP is positive. With respect to iPohp{o'q) 
in simulations and HNC/MSA, from Fig. [5] it is evident 
that, for those 1:1 systems, the value (and sign !) of the 
"^o/zp-integral is determined by overcharging, for (Tq > 0, 
and by charge reversal, for ao < 0. For positive colloidal 
charges, the increment of the height of the overcharging 
peak at OHP when zm augments dictates the increas- 
ing positiveness of ^ohp- In addition, we notice that, 
for negative macroparticle charges, a dominant peak of 
charge reversal occurs either at the OHP or beyond that 
point, which in the present work is a distinctive feature 
of those descriptions that incorporate consistently the 
ionic size and size asymmetry effects (i.e. simulations 
and HNC/MSA) and, therefore, it is absent in URMGC. 
Such strong charge reversal gives rise to the positive sign 
of ipOHP for negatively charged colloids, and then com- 
pletes our explanation of the fact that, for simulations 
and HNC/MSA, ipoHP is always positive for the range 
covered in Fig. [7)3. 

In connection with electrophoresis experiments, and 
depending on the location of the shear plain, our results 
suggest two possible scenarios. If the sign of the elec- 
trophoretic mobility, ^, were associated to the sign of 
the potential ipiHP, our treatment of the size- asymmetric 
EDL would predict the inversion of the colloidal mobil- 
ity very near to the point of zero charge. In contrast, if 
the sign of /x were associated to the sign of the potential 
ipOHP-, our study would foresee positive mobilities in a 
wider range of colloidal charges around zero. In addi- 
tion. Fig. [7)3 resembles a situation reported in Fig. 2b 



of Ref. [231 by Johnson et al., where zeta potential mea- 
surements of a-alumina in presence of 1 M LiNOa are 
plotted as a function of pH. In those results, the zeta 
potential is always positive in a wide interval of pH that 
encompasses the point of zero charge, which parallels our 
results in Fig. [7)). This would indicate that such positive 
values of the zeta potential are due in part by the pref- 
erential adsorption of one of the species (Li+), which in 
turn should be induced by the ionic size asymmetry, as 
deduced from our survey of the spherical EDL. Although 
other complex mechanisms are into play in these kind of 
experiments [sl, HO, [zl, [ill , our size-asymmetric model 
seems to capture adequately relevant phenomena occur- 
ring in real systems, and hence it could be considered as 
a basic representation of the EDL to which further im- 
provements (as Van der Waals dispersion forces or more 
sophisticated chemical mechanisms) can be incorporated 
in order to predict experimental data more accurately. 



B. Divalent size-asymmetric electrolytes 

We proceed to investigate the properties of the electri- 
cal double layer for a macroion immersed in a 2:2 salt, i.e. 
for systems with stronger charge correlations. As it will 
be shown, many of the findings reviewed in the prior sec- 
tion devoted to univalent electrolytes are also present, in 
an enhanced way, for the case of divalent ions. As before, 
we examine first the EDL at the point of zero charge and 
later we ponder instances with charged macroparticles. 
In what follows, we consider EDL systems formed by a 
colloid and a bath of a 2:2, 0.5 M electrolyte, with the 
same diameters as specified in the earlier section. 

The simulational, HNC/MSA and URMGC radial dis- 
tributions of divalent ions around a macroparticle with 
Zm = are included in Fig. [8| (main panel). There 
we have incorporated the Monte Carlo and HNC/MSA 
pair correlation functions for the associated hard-sphere 
mixture (uncharged "ions"). Again, the idea is to exem- 
plify how the coulombic correlations modify the struc- 
ture of the pure hard-sphere fluid, in order to gain some 
insight into the relative importance of the entropic and 
charge correlations. From the direct contrast between 
the simulational RDFs of hard spheres and ions in the 
main panel of the figure, we observe that the charge ef- 
fects are very strong, completely modifying the accumu- 
lation of the smaller species around the macroparticle. 
In fact, the changes in the Monte Carlo gj{r') of divalent 
ions with respect to the RDFs of hard spheres are bigger 
than those occurring in univalent systems (compare the 
main panels of Figs. [T]and[8|). For the simulations of 
2:2 electrolytes, the impact of the valence is much more 
important for the smaller ionic species, where we observe 
that the corresponding contact population depletes so 
much that the concentration of cations near the surface 
is below the bulk value. A similar dewetting has been 
reported in a recent study of the planar EDL for size- 
asymmetric 1:1, 1:2 and 1:3 salts |53| , with the multiva- 
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lent ions corresponding to the smaller species. In that 
work it was found that the amount of small cations in 
the proximities of a neutral plane decreased as long as 
the electrostatic coupling (i.e. the valence) was height- 
ened, being our results consistent with such behavior. On 
the other hand, in Fig. [8] (main panel) we can appreci- 
ate that the HNC /MSA results for hard spheres compare 
well with the simulations, whereas for the divalent case 
this formalism overestimates the Monte Carlo data for 
the two ionic species. Furthermore, HNC/MSA fails to 
describe the correct tendency for the RDF of big anions, 
predicting an increase of such function with respect to 
the corresponding RDF of the larger hard spheres, which 
is clearly not the behavior seen in the MC data. De- 
spite this poor achievement in the estimation of the gj{r') 
for the electrolyte, HNC/MSA does well at the level of 
the local electrolyte charge density (see inset of Fig. [8|), 
capturing the correct trends of the simulations. This is 
not a surprising circumstance given that the dependence 
of p*{r') on the difference g+{r') — g~{r') explains why 
the notable discrepancies existing between the RDFs of 
Monte Carlo and HNC/MSA are attenuated when the 
local charge densities are compared instead. URMGC, 
on the other hand, highly disagrees with the simulations 
for the gj{r') and p*{r'). Special note must be taken 
of the very low URMGC concentration of small ions in 
the Helmholtz zone, which, analogously to the univalent 
case, will determine the weak intensity of the CR and 
OC phenomena for 2:2 systems in this theory. 

The simulational and theoretical ionic distributions of 
the divalent salt yield the accumulated charge and mean 
electrostatic potential given in Fig. [5] and its inset. As 
it was detected in the monovalent situation, the ionic 
size asymmetry promotes an adsorbed layer of charge 
and the concomitant existence of a zero surface charge 
double layer (see Fig. [5]), as well as a non-zero MEP 
at (and between) the Helmholtz planes. The adsorbed 
charges, up to the OHP, are w 2.6e, 3.3e, and 2.0e for 
MC, HNC/MSA, and URMGC, respectively, which are 
smaller than those observed for the 1:1 salt. In turn, con- 
trasted with Monte Carlo, the MEP within the Helmholtz 
planes is overestimated by HNC/MSA and underesti- 
mated by URMGC. Globally we see that the performance 
of HNC/MSA and URMGC is very similar to that dis- 
played for the monovalent salt, with a better qualitative 
similitude between HNC/MSA and the simulations. For 
URMGC the P{r') and '0(r') profiles evidence the re- 
striction of all the size asymmetry and hard-core effects 
to happen uniquely in the Helmholtz zone. 

The Monte Carlo and theoretical RDFs, LECD, IC, 
and MEP curves corresponding to the divalent system 
and two values of the colloidal charge, namely zm — —4 
and Zm = 4, are plotted in Figs. [TU] and [TT] When 
Zm = —4 (Fig. [TUb). there is a strong adsorption of coun- 
terions inside the Helmholtz planes, accompanied by an 
important depletion of colons, particularly for MC and 
HNC/MSA (main panel). Besides, it is worth noticing 
that very close to the OHP the concentration of colons is 



smaller than that of counterions (compare with Fig. [3^), 
which seems to indicate that at these distances the repul- 
sion forces between the macroion and the colons are far 
from being screened. The associated local charge densi- 
ties are portrayed in the inset, where we have curves that 
do not jump to a negative value immediately after the 
OHP, showing that the macroion-coions repulsion propa- 
gates beyond the Helmholtz planes. When zm = 4 (Fig. 
ITOb). and despite of the intense electrostatic repulsion, 
the MC, HNC/MSA and URMGC concentrations of the 
small colons close to the IHP turn out to be small but 
finite, effectively increasing the charge of the macroion, 
and therefore driving a significative enhancement in the 
counterions concentration near the OHP (main panel). 
As a result, the LECD profiles develop a small first pos- 
itive layer followed by an ample negative region (see in- 
set). 

In Fig. [TT]we present the integrated charges and MEPs 
for the divalent salt and macroion valences zm = —4 
and Zm = 4. For zm = —4 (panel (a)), the ICs of 
Monte Carlo and HNC/MSA exhibit charge reversal in 
the Helmholtz zone, but the maximum charge reversals 
are outside that region (at r' /D+ « 1.5 and r' /D+ w 1.4 
for simulations and HNC/MSA, respectively). The cor- 
responding MEPs, on the other hand, display a strong 
reversal, with the maximum potential being inside the 
Helmholtz planes (see inset). In the case zm = 4 (panel 
(b)) we see that the MC and HNC/MSA integrated 
charge profiles present overcharging, followed by alter- 
nate oscillations of decreasing amplitude, whereas the 
associated MEPs (in the inset) decay in a fluctuating 
manner to zero. This observation of OC in the computer 
"experiments" of divalent systems complements our pre- 
vious findings for univalent salts and, thus, consolidates 
our simulational proof of overcharging as a genuine fea- 
ture of size- asymmetric primitive model EDLs. 

A simultaneous analysis of all the structural infor- 
mation for divalent systems contained in Figs. [TO] and 
[Illshow that the HNC/MSA theory follows closely the 
Monte Carlo data, whereas the URMGC approach ex- 
hibits notable differences with respect to the simulations. 
From these figures it is also verified that in the URMGC 
description of the EDL there is a confinement of the 
steric-related phenomenology (e.g. CR, OC and the non- 
monotonic character of the structural functions) within 
the Helmholtz zone. 

Fig. [12] portrays the Monte Carlo and HNC/MSA in- 
tegrated charges for varying zm (—8 < zm < 8) in or- 
der to analyze the evolution of the size asymmetry ef- 
fects when the divalent systems depart from the point 
of zero colloidal charge. For zm < —2 (Fig. \12k ) the 
IC profiles of MC and HNC/MSA present charge rever- 
sal, with the maximum located outside the Helmholtz 
zone. As the charge increases (towards zm — 0) the 
location of the peak of each curve shifts to shorter dis- 
tances and eventually lands at the OHP. Further increase 
of Zm leads to the results plotted in Fig. [T^ . where the 
system displays overcharging, with a maximum intensity 



(i.e. \P{r') — zm\) that decreases with the augment of 
the macroion's valence. This overcharging phenomenon 
only takes place in between the Helmholtz planes, and 
it is followed by a first region of charge reversal. At the 
limiting valence zm ~ 16 overcharging virtually disap- 
pears; meanwhile the charge reversal beyond the OHP 
continues to exist with an increasing magnitude. All 
these characteristics of P{r') for zm > are equally de- 
scribed for both the simulations and HNC/MSA. The 
complete development of the P(r') for 2:2 electrolytes 
and — 8 < zm < 8 reveals a passage from charge rever- 
sal to overcharging. This crossover at zm = is due 
uniquely to the size asymmetry of the salt ions, as the 
change of sign of zm merely inverts the role of the anions 
and cations (as colons or counterions). 

Therefore, in the present investigation it is has been 
found that the charge reversal and overcharging observed 
near the point of zero charge are mainly caused by en- 
tropic contributions coming from the size-asymmetric 
nature of the electrolyte ions. However, when zm is 
increased, the charge reversal and overcharging in the 
Helmholtz zone disappear, and only CR persists. This 
remaining charge reversal beyond the OHP comes from 
the interplay between excluded volume and electrostatic 
correlations and, thus, it is not exclusive of a size- 
asymmetric model. That is the reason why charge re- 
versal has been already reported in many studies of the 
restricted primitive model EDL, in which the excluded 
volume effects are also consistently taken into account 
[T^ I . Accordingly, EDL theories for genuine punctual ions 
and uniformly charged surfaces do not predict charge re- 
versal at all. 

We would like to mention that our divalent systems 
also exhibit regions of local inversion of the derivatives 
of the integrated charge (LIDIC) . These zones are delim- 
ited in Figs. [T^ andlT^ by the upward arrows for Monte 
Carlo and by the downward arrows for HNC/MSA. In- 
side these "anomalous" regions the accumulated charge 
at any point decreases when the colloidal charge is aug- 
mented and vice versa. Similarly to the univalent case, 
and starting from the definition of the MEP (see Eq. ([8])), 
the existence of this LIDIC domains provides a way to 
explain the behavior of the potential-charge relationship 
at the Helmholtz planes. 

The corresponding sequence (—8 < zm < 8) of UR- 
MGC integrated charges for 2:2 electrolytes is given in 
Fig. [T31 Visibly, each one of the profiles in this series 
variates monotonically with respect to the distance r' 
and, as a whole, the full set of curves changes monotoni- 
cally as a function of zm- In other words, after the OHP, 
every P(r') goes to zero without spatial oscillations, and 
the complete sequence does not show LIDIC regions. 

The MEP of MC, HNC/MSA, and URMGC at the 
Helmholtz planes are shown in Fig. [T3] as a function of 
the surface charge ctq- At the IHP, all the approaches 
exhibit a MEP that is monotonic and displays a posi- 
tive zero-charge potential and a negative ctq''** such that 
ipiHP > for erg'"** < (To < 0, just hke in the monova- 
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lent systems. From panel (b) we also find that tpOHP is 
also monotonic for URMGC. In contrast to the mono- 
valent case, however, ipoHP from MC and HNC/MSA 
is non-monotonic and displays a minimum at a positive 
(T™"\ This augment of tpoHP when ctq decreases below 
fjmm gygggsts that there could be experimental systems 
in which the reverted electrophoretic mobility always in- 
creases its magnitude. 

The trends of these MEPs can be explained in terms 
of the corresponding IC profiles (see Figs. [T^ and [T5)) . as 
discussed previously in the context of monovalent elec- 
trolytes. In particular, the concavity of ipoHP for MC 
and HNC/MSA can be explained by observing in Fig. 
[12b the behavior of the integral of P(r)/r^ between two 
consecutive electroneutral points of P(r) , and noting that 
the absolute value of these negative areas decreases when 
zjii diminishes. Thus, for ctq > 0, the area from the OHP 
to the first electroneutral point is positive and larger than 
the negative area between the first and second electroneu- 
tral points (see Fig. WBd). Since the difference between 
these two contributions to the ipoHP increases with ctq, 
this explains the positive value of ipoHP and its increas- 
ing behavior. An analogous argument explains the trends 
observed when ao < (Fig. [T^). despite the disappear- 
ance of the maximum charge reversal at the OHP when 
Zm decreases. 

The concavity of 4'Ohp can also be cast in terms of 
the local inversion of derivative of the MEP (LIDMEP), 
which we define analogously to the LIDIC. With this in 
mind. Fig. [TSl shows the MEP of the divalent system for 
several values of zm- For zm < the MC data shows 
a region of LIDMEP in 0.82 < r'/D+ < 2.6 (upward 
arrows in Fig. I15b .l. A similar region is seen in Fig. 115b 
(zm > 0). The downward arrows in these figures indicate 
the location of such regions for HNC/MSA. Therefore, 
the change in sign of the derivative of the MEP observed 
in Fig. [Hb appears because the OHP is in a LIDMEP 
region, where 



^^(""'J^ = < 0. (13) 

Since ionic size and size asymmetry are evidently im- 
portant in real highly coupled electrokinetic systems, 
these last results suggest that the interpretation of the 
zcta potential in electrokinetic experiments must be 
done carefully because the experimental MEP profiles 
could present LIDMEP regions, according to our model 
calculations. Consequently, the measured C might corre- 
spond to potentials located in LIDMEP regions, giving 
rise to the possibility of decreasing C for increasing (Tq 
and vice versa. The occurrence of these non-monotonic 
effects in the potential are clearly precluded in the PB- 
based interpretation of the experimental measurements, 
as evidenced in the present study. 
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V. SUMMARY AND CONCLUSIONS 

In this paper the size-asymmetric electrical double 
layer (EDL) of 1:1 and 2:2 salts around a slightly charged 
spherical macroparticle was studied by assuming the 
primitive model of an electrolyte as a representation of 
the ionic bath, and using Monte Carlo simulations and 
the HNC/MSA integral equation in order to calculate the 
corresponding properties of the model system. Therefore, 
in the simulation and HNC/MSA integral equation the 
important ionic size and size asymmetry effects have been 
fully incorporated through the consistent consideration 
of the colloid-ions and ion-ion interactions. Additionally, 
the possible consequences of a partial treatment of such 
excluded volume and size asymmetry contributions in the 
properties of the spherical EDL were assessed via the 
comparison with the predictions of the classic URMGC 
theory for spherical geometry. In the semi-punctual UR- 
MGC formalism the finite nature of the ions is considered 
only through the use of two different distances of closest 
approach between the colloid and counterions and col- 
loid and colons. Our Monte Carlo data evince that the 
finite size and size asymmetry of the ionic species near 
a barely charged colloid produce remarkable phenomena, 
such as charge reversal, overcharging, the existence of 
potentials of zero charge, and of a zero charge electri- 
cal double layer, as well as the possibility of oscillations 
in the ionic densities, and of sign alternacy in the mean 
electrostatic potential. The simulations also show that 
the extent of all these features is intimately linked to the 
ability of the smallest ions to penetrate the Helmholtz 
zone. Importantly, in this report we present the first 
simulational corroboration of overcharging (i.e. of the 
increment of the native colloidal charge prompted by an 
anomalous adsorption of colons), a peculiarity originally 
predicted by Jimenez- Angeles et al. (23| in a more com- 
plex EDL system. Here it is evidenced that, close to 
the point of zero charge, overcharging is a direct con- 
sequence of the size asymmetry of the ions. On the 
whole, the HNC/MSA results coincide with the simula- 
tions, then providing an essentially correct picture of the 
size- asymmetric SEDL near the point of zero charge. In 
contrast, URMGC disagrees quantitatively and qualita- 
tively with the Monte Carlo trends in most of the situa- 
tions examined here. Notably, and due to its inconsistent 
treatment of the hard-core and electrostatic correlations, 
URMGC exhibits a spatial localization of all its ionic 
size and size asymmetry effects within the zone delim- 
ited by the inner and outer Helmholtz planes. One of 
the most characteristic sequels of the full incorporation 



of the ionic correlations in all the extension of the EDL, 
and not only in the Helmholtz zone, is the occurrence 
of a non-monotonic behavior in the ionic concentration, 
charge density and mean electrostatic potential profiles 
of Monte Carlo and HNC/MSA. In this respect, the ex- 
istence of spatial regions in which the integrated charge 
and mean electrostatic potential present a non-uniform 
comportment with respect to the variation of the colloidal 
charge proves its relevance for the understanding of the 
potential-charge relationship at the Helmholtz planes. In 
particular, such "non-uniform" regions serve to explain 
the existence of a singular minimum in the ipoHP — ctq 
curve. 

The plausible identification between the well-known 
zeta potential and the mean electrostatic potential 
around the Helmholtz zone, which is a usual hypothesis 
in the interpretation of electrophoretic mobility measure- 
ments, leads us to suggest several phenomena (e.g. the 
presence of an anomalous sign in the zeta potential and 
the occurrence of increased reversed mobilities) which 
could be the objectives of future experimental protocols. 
As a first attempt in this direction, we have pointed out 
the consistency between our results and some experimen- 
tal data of electrokinetic mobilities for a-alumina parti- 
cles. The eventual confirmation in the laboratory of any 
of the theoretical predictions presented in this work (for 
example, charge reversal and overcharging) by means of 
some electrokinetic or static technique could indicate the 
pertinence of our simple primitive-model-based represen- 
tation of the size-asymmetric spherical EDL as a start- 
ing point to develop more faithful descriptions of real 
colloidal systems, either in equilibrium or under exter- 
nal fields. Finally, and in relation with this, the rich 
phenomenology discussed in the present communication, 
that arises from both the ionic size and size asymmetry 
effects, puts a word of caution about the possible usage 
of the Poisson-Boltzmann viewpoint and other contingent 
formalisms, such as the so-called standard electrokinetic 
model [t^ , outside their range of applicability. 
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FIG. 1: Radial distribution functions of a size-asymmetric bi- 
nary mixture of hard spheres and a size-asymmetric 1:1 elec- 
trolyte, both at a 1 M concentration, around an uncharged 
{zm = 0) macroparticle (in the main panel). Thereinafter, the 
diameters of the ionic species and colloid are D+ = 4.25 A, 
D- = 8.5 A, and Dm = 20 A, respectively. The small 
hard spheres have the same diameter as that of the cations, 
D\ = 75+ , and the large hard spheres diameter is the same as 
that of the anions, D2 = D- . The open circles correspond to 
Monte Carlo (MC) simulations of the size-asymmetric mix- 
ture of hard spheres and the open squares represent the MC 
data of the size-asymmetric monovalent electrolyte. The solid 
lines are associated to the asymmetric mixture of hard spheres 
in the HNC/MSA approach. The dashed and dot-dashed 
lines are the HNC/MSA and URMGC results for the 1:1 
size-asymmetric salt, respectively. In the inset we depict the 
corresponding local electrolyte charge density for the 1:1 elec- 
trolyte case. Open triangles, and dashed and dot-dashed lines 
correspond to simulations, and the HNC/MSA and URMGC 
theories, respectively. Here and in the rest of the figures, the 
distance r' is measured from the colloidal surface. 
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FIG. 2: Integrated charge (main panel) and mean electro- 
static potential (inset) of a size-asymmetric 1:1, 1 M elec- 
trolyte around an uncharged {zm ~ 0) macroparticle. The 
open circles, solid and dot-dashed lines correspond to the MC, 
HNC/MSA and URMGC data, respectively. 
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FIG. 3: Radial distribution functions of a size-asymmetric 1:1, 
1 M electrolyte, around a charged macroion (main panel). 
In Fig. EJa) zm = -4, and in Fig. E^b) zm = 4. The 
open squares, solid and dot-dashed lines represent the MC, 
HNC/MSA and URMGC results, respectively. In the insets 
the corresponding local electrolyte charge densities are por- 
trayed. Open triangles, and solid and dot-dashed lines are 
associated to simulations, and the HNC/MSA and URMGC 
theories, respectively. 
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FIG. 4: Integrated charge (main panel) and mean electro- 
static potential (inset) of a size-asymmetric 1:1, 1 M elec- 
trolyte around a charged macroion. In Fig. |4ja) zm ~ —4, 
and in Fig. IHb) zm ~ 4. The open circles, solid and dot- 
dashed lines correspond to the MC, HNC/MSA and URMGC 
data, respectively. 
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FIG. 5: Integrated charge of a size- asymmetric 1:1, 1 M elec- 
trolyte around a charged macroion. In Fig. Efa) the open 
circles, squares, diamonds and triangles correspond to Monte 
Carlo simulations for zm ~ —2, —4, —6, —8, respectively, and 
in Fig. [5jb) are associated to zm = 2,4,6,8, respectively. 
In Fig. [5ja) the solid, dashed, dot-dashed and dotted lines 
correspond to HNC/MSA resuhs for zm = -2, -4, -6, -8, re- 
spectively, and in Fig. [5jb) are associated to zm ~ 2,4,6,8, 
respectively. 
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FIG. 6: The same as in Fig. [5] but for the URMGC theory. 
In the main figure the solid, dashed, dot-dashed and dotted 
lines correspond to zm = —2, —4, —6, —8, respectively, and in 
the inset to zm ~ 2, 4, 6, 8, respectively. 




FIG. 7: Mean electrostatic potential (MEP) at the Helmholtz 
planes as functions of the surface charge density of a macroion 
immersed in a size-asymmetric 1:1, 1 M electrolyte. The open 
circles, and solid and dashed lines correspond to simulations, 
and the HNC/MSA and URMGC theories, respectively. In 
Fig. [3; a) the MEP at the IHP is plotted, and in the Fig. [T^b) 
the same is done for the MEP at the OHP. 
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FIG. 8: Radial distribution functions of a size-asymmetric 
binary mixture of hard spheres and a size-asymmetric 2:2 
electrolyte, both at a 0.5 M concentration, around an un- 
charged (zm = 0) macroparticle (in the main panel). The 
small hard spheres have the same diameter as that of the 
cations, Di = D+, and the large hard spheres diameter is 
the same as that of the anions, D2 = -D-. The open cir- 
cles correspond to Monte Carlo (MC) simulations of the size- 
asymmetric mixture of hard spheres and the open squares 
represent the MC data of the size-asymmetric divalent elec- 
trolyte. The solid lines are associated to the asymmetric mix- 
ture of hard spheres in the HNC/MSA approach. The dashed 
and dot-dashed lines are the HNC/MSA and URMGC results 
for the 2:2 size-asymmetric salt, respectively. In the inset we 
depict the corresponding local electrolyte charge density for 
the 2:2 electrolyte case. Open triangles, and dashed and dot- 
dashed lines correspond to simulations, and the HNC/MSA 
and URMGC theories, respectively. 
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FIG. 9: Integrated charge (main panel) and mean electro- 
static potential (inset) of a size-asymmetric 2:2, 0.5 M elec- 
trolyte around an uncharged {zm = 0) macroparticle. The 
open circles, solid and dot-dashed lines are associated to the 
MC, HNC/MSA and URMGC data, respectively. 
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FIG. 10: Radial distribution functions of a size-asymmetric 
2:2, 0.5 M electrolyte, around a charged macroion (main 
panel). In Fig. [lOja) zm = -4, and in Fig. ^[h) zm = 4. 
The open squares, solid and dot-dashed lines represent the 
MC, HNC/MSA and URMGC resuhs, respectively. In the 
insets the corresponding local electrolyte charge densities are 
portrayed. Open triangles, and solid and dot-dashed lines are 
associated to simulations, and the HNC/MSA and URMGC 
theories, respectively. 
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FIG. 11: Integrated charge (main panel) and mean electro- 
static potential (inset) of a size-asymmetric 2:2, 0.5 M elec- 
trolyte around a charged macroion. In Fig. IllT a) zm = —4, 
and in Fig. Illf b') zm ~ 4. The open circles, solid and dot- 
dashed lines correspond to the MC, HNC/MSA and URMGC 
data, respectively. 
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FIG. 12: Integrated charge of a size-asymmetric 2:2, 0.5 M 
electrolyte around a charged macroion. In Fig. I12f a') the open 
circles, squares, diamonds and triangles correspond to Monte 
Carlo simulations for zm ~ —2, —4, —6, —8, respectively, and 
in Fig. I12f b) are associated to zm ~ 2, 4, 6, 8, respectively. 
In Fig. I12f a) the solid, dashed, dot-dashed and dotted lines 
correspond to HNC/MSA results for zm = —2, —4, -6, -8, re- 
spectively, and in Fig. I12f bl are associated to zm ~ 2, 4, 6, 8, 
respectively. 
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FIG. 13: The same as in Fig. [T2]but for the URMGC theory. 
In the main figure the solid, dashed, dot-dashed and dotted 
lines correspond to zm = —2, —4, —6, —8, respectively, and in 
the inset to zm ~ 2, 4, 6, 8, respectively. 




FIG. 14: Mean electrostatic potential (MEP) at the 
Helmholtz planes as functions of the surface charge density 
of a macroion immersed in a size-asymmetric 2:2, 0.5 M elec- 
trolyte. The open circles, and solid and dashed lines corre- 
spond to simulations, and the HNC/MSA and URMGC the- 
ories, respectively. In Fig. [HJa) the MEP at the IHP is 
plotted, and in the Fig. I14f bl the same is done for the MEP 
at the OHP. 




FIG. 15: Mean electrostatic potential ol a size-asymmetric 
2:2, 0.5 M electrolyte around a charged macroion. In Fig. 
I15r a) the open circles, squares, diamonds and triangles corre- 
spond to Monte Carlo simulations for zm ~ —2, —4, —6, —8, 
respectively, and in Fig. I15f b) are associated to zm = 
2, 4, 6, 8, respectively. In Fig. I15f a) the solid, dashed, dot- 
dashed and dotted lines correspond to HNC/MSA results for 
Zm = —2, —4, —6, —8, respectively, and in Fig. I15r b) are as- 
sociated to Zm = 2, 4, 6, 8, respectively. 



